Low-Energy Inelastic Neutrino Reactions on 4 He 
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The inelastic scattering of neutrino off 4 He is calculated microscopically at energies typical for 
core collapse supernova environment. The calculation is carried out with the Argonne V18 nucleon- 
nucleon potential and the Urbana IX three nucleon force. Full final state interaction is included via 
the Lorentz integral transform (LIT) method. The contribution of axial meson exchange currents 
to the cross sections is taken into account from effective field theory of nucleons and pions to order 
0(Q 3 ). 
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The current theory of core collapse supernova holds 
some open questions regarding the explosion mechanism 
and late stage nucleosynthesis. In order to analyze these 
questions, a better understanding of the involved micro- 
scopical processes is needed. In particular, due to the 
high abundance of a particles in the supernova environ- 
ment, the inelastic neutrino- 4 He reaction has drawn at- 
tention in recent years. This interest yielded a number of 
studies trying to estimate the cross-section and the role 
of neutrino- 4 He reactions in the described phenomena, 
[E HIE H H, IE 0, S|- However to date, a full ab-initio 
calculation that includes a realistic nuclear Hamiltonian 
is still missing. Moreover, the contribution of meson ex- 
change currents (MEC) to this particular scattering pro- 
cess was never estimated. 

In this letter we present a full ab-initio calculation 
of the inelastic neutrino- 4 He reactions that meets these 
challenges. Specifically, we consider the energy de- 
pendent inclusive inelastic cross-sections for the follow- 
ing channels ^e^,^)^, 4 rIe(j> a ;,j> :E ')|X, 4 He(i/ e ,e + ) 4 X, 
and 4 He(^ e ,e~)|X, where x = e,/i,r and stands for 
the final state ^4-nucleon system, with charge Z. 

Core collapse supernovae are believed to be neutrino 
driven explosion of massive stars. As the iron core of the 
star becomes gravitationally unstable it collapses until 
the nuclear forces halt the collapse and drive an outgoing 
shock. This shock gradually stalls due to energy loss 
through neutrino radiation and dissociation of the iron 
nuclei into a mixture of a particles and free nucleons. 

At this stage, the proto-neutron star cools mainly by 
emitting neutrinos in enormous numbers. These neutri- 
nos are a result of thermal pair production, and thus are 
produced in flavor equilibrium. The characteristic tem- 
peratures of the emitted neutrinos are about 6 — 10 MeV 
for v a ,T (*V,t)j 5—8 MeV for D e , and 3—5 MeV for v e . The 
difference in temperature originates from the large cross- 
sections for v e , v e electron scattering and charge current 
reactions. 

At these temperatures there is a considerable amount 
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of [i and r neutrinos (and anti-neutrinos) which carry 
more than 20 MeV. At such energies the v^r's (P^r's) 
can dissociate the 4 He nucleus through inelastic neutral 
current reactions, thus creating the seed to light element 
nucleosynthesis in the supernova environment. A knock 
out of a nucleon from a 4 He nucleus in the helium rich 
layer, followed by a fusion of the remaining trinucleus 
with another a particle, will result in a 7-body nucleus. 
This process is an important source of 7 Li, and of 11 B 
and 19 F through additional a capture reactions. Due to 
the high dissociation energy of the a, this mechanism is 
sensitive to the high-energy tail of the neutrinos. Thus a 
correct description of the process must contain an exact, 
energy dependent cross-section for the neutral inelastic 
a — v reaction, which initiates the process. 

The relatively low temperature of the z/ e 's and D e 's 
emitted from the star's core suppress the probability for 
inelastic reactions of these neutrinos with 4 He in the su- 
pernova scenario. Oscillations of the \x and r (anti) neu- 
trinos can yield a secondary source of energetic electron 
neutrinos. Apparently, there is a 1-3 neutrino oscillation 
resonance in O/C layer. As a result, substantial incre- 
ment in the charged current reaction rates, and thus in 
the abundance of 7 Li and 11 B, can take place in the he- 
lium layer, depending on the value of the mixing angle 
013 . In fact, such an effect can be used to constrain the 
value of this oscillation parameter Q • 

The possible role of inelastic v — ot reactions in reviving 
the supernova explosion shock was pointed out by Haxton 
[l[ . The hot dilute gas above the proto-neutron star and 
below the accretion shock contains up to 70% 4 He nuclei. 
It is believed that neutrinos emitted from the collapsed 
core (the proto-neutron star) deposit energy in the mat- 
ter behind the shock, and eventually reverse the flow and 
revive the shock. This delayed shock mechanism, origi- 
nally introduced by Colgate and White Q, has not yet 
been proved in full hydro-reactive simulations. Haxton 
has suggested that inelastic neutral reactions of neutri- 
nos with 4 He can lead to an enhanced neutrino energy 
deposition. This effect is usually ignored (see however 
[E [8j] ) and was not yet considered in a full hydrodynamic 
simulation. 

The first challenge in the study of the inelastic 
neutrino- 4 He reactions is the solution of the four-body 
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problem, for ground and excited states. As the 4 He has 
no bound excited states, a detailed knowledge of the four 
nucleons continuum is needed to assure the final state 
interaction (FSI). This makes an explicit calculation im- 
possible, since a complete description of the nuclear four- 
body system is currently out of reach. We avoid this 
complication by calculating the FSI through the Lorentz 
integral transform (LIT) method For the solution 

of the ground state wave function and the LIT equa- 
tions we use the effective interaction in the hyperspher- 
ical harmonics (EIHH) approach [ll|, [12, EH- For the 
nuclear Hamiltonian, we take the nucleon-nucleon (NN) 
potential Argonne V18 (AV18) [lj] with the Urbana IX 
(UIX) [H[ three nucleon force (3NF). This Hamiltonian 
has been used successfully to reproduce the spectra of 
light nucl ei [ Lai, and electro-weak reactions with light 
nuclei p^.ll7l.ll8j. 

In the limit of small momentum transfer with respect 
to the mass of the Z 1 W ± bosons, the weak interaction 
Hamiltonian is given by H\y — — ^ / d^xj^x^jJ^ix), 

where G is the Fermi weak coupling constant, j^ix) is the 
lcpton current, and is the nuclear current. The lepton 
is a point Dirac particle, and evaluating its current and 
its contribution to the cross-section is relatively simple, 
yielding only kinematical factors. The nuclear current, 
however, is more complicated. The formal structure of 
the nuclear weak neutral current is, 

J° = (l-2.sin 2 My^ + y^-2-sin 2 ^ij / r , (1) 
and the structure of the charged currrents is, 

% = f # + T J t ■ W 

where the superscript A (V) stands for axial (vector) 
currents. 

The leading contributions to these operators are the 
one-body terms. It is well known, however, that mesonic 
degrees of freedom can contribute to the nuclear currents 
through many-body terms, namely MEC, even if they 
do not appear explicitly in the Hamiltonian. The mod- 
ern point of view [HI has created a systematic way of 
considering these degrees of freedom, that is the effec- 
tive field theory (EFT) approach. EFT is based on the 
idea that an observable characterized by a momentum Q, 
does not depend on momenta much higher than Q. One 
introduces a cutoff momentum A, and integrates out the 
degrees of freedom present at Q larger than A. A per- 
turbation theory in the small parameter Q/A can now 
be developed systematically. The coefficients of the dif- 
ferent terms are called low-energy constants (LEC), usu- 
ally calibrated in experiments or by theory. EFT has two 
major advantages, one is the link to the underlying high- 
energy theory, which in the case of the strong interaction 
is commonly believed to be Quantum Chromodynamics 
(QCD). The other advantage is the ability to provide a 
control of the accuracy in the calculation. The problem 



with EFT is that while a percentage level accuracy in 
describing scattering process is already achieved using a 
next-to-leading order (NLO) lagrangian, this is not the 
case when trying to recover successfully the wealth of ex- 
perimental data described by the phenomenological ap- 
proach (nuclear binding energies for example) . This task 
demands at least next-to-next -to-next-to-leading order 
(N 3 LO) EFT lagrangians @, HI- It is thus clear that 
a hybrid method that joins together the success of the 
standard nuclear physics approach and the clear advan- 
tages of EFT is called for, although the resulting MEC 
will not be completely consistent with the nuclear Hamil- 
tonian. This hybrid approach was coined by Rho as 
MEEFT ("more effective EFT") [H, and was already 
applied to study electro- weak reactions for A = 2,3,4 
nuclei 23]. In this work we adopt the hybrid approach 
combining the phenomenological AV18 and UIX nuclear 
potentials with EFT based nuclear MEC. 

The conservation of vector current (CVC) hypothesis 
states that the vector current is an isospin rotation of 
the electromagnetic current. Thus, the electric part of 
the vector meson exchange currents can be approximated 
very well at low q via the Siegert theorem, from the single 
nucleon vector charge operator. That is not the case for 
the axial current, which is not conserved and should be 
calculated explicitly. For this task we shall use the EFT 
meson exchange currents. 

The typical energy scale of the neutrino in the super- 
nova environment is some tens of MeV, thus a proper 
cutoff is of at least a few hundred MeV. It is important to 
notice that the cutoff should not be higher than the mass 
of the nucleon, which is the order of the QCD mass scale. 
We will use cutoff values in the range A = 400 — 800 MeV. 
In the EFT scheme employed here, nucleons and pions 
are the explicit degrees of freedom. The model includes 
the pions as goldstone bosons of the chiral symmetry [24[ • 
The axial currents are the Nother currents derived from a 
NLO order lagrangian, in a relativistic approach. These 
currents are accurate to N 3 LO and are given in momen- 
tum space, as they originate from a Lorentz invariant 
theory. For the transformation to configuration space, 
we perform a Fourier transform with cutoff [23| . 

6{x) = J -0^6(q)S A (q). (3) 

The cutoff function Sa(q) is 1 for q <C A, and approaches 
for g ^ A. We use a gaussian cutoff function as pro- 
posed by Park et. al. [23|, i.e. Sa(<z) = exp(-jj). It 
is important to note that this method leads to the same 
single nucleon operators as the standard nuclear physics 
approach. The meson exchange currents in configura- 
tion space are the Fourier transform of propagators with 
a cutoff, which in the limit A — > oo are just the usual 
Yukawa functions. In contrast to the standard nuclear 
physics approach, the coefficients of the functions are not 
structure functions, but LECs. All the LECs which orig- 
inate from a nucleon-pion interaction, namely f v , 
C3, and C4, are calibrated using low-energy pion-nucleon 
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scattering. Alas, in this order there are also two nucleon 
contact terms, which introduce LECs that can be cali- 
brated only by nuclear matter processes. Fortunately, us- 
ing Lorentz invariance, the axial currents introduce only 
one unknown LEC, which is denoted by d r . This coef- 
ficient can be calibrated by the triton half life time for 
each value of A, thus d r — d r (A). In our calculations we 
reproduce the results published by Park et. al. [23j], for 
the cutoff dependency. 

For low-energy reactions a multipole decomposition of 
the currents is useful. Applying Fermi's golden rule, to 
inclusive reactions with unpolarized targets, and consid- 
ering recoil effects, the differential cross-section takes the 
form [U, 
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where kf is the momentum of the outgoing lepton, J L = 
is the angular momentum of the 4 He, Zf is the charge of 
the residual nuclear system. The four- vector (to, q) repre- 
sents energy and momentum transfer, and is the angle 
between the incoming neutrino direction and outgoing 
lepton direction. The superscript a denotes the isospin 
component, with a = for the neutral current and a = ± 
for the charged currents. The Coulomb factor F a (Z,k) 
is equal to 1 for neutral currents, and is the Fermi func- 
tion for charged current. The functions Xdi6 2 are the 
leptonic kincmatical factors (related to the 0\ , O2 multi- 
poles, Xq i = Xq i q i ). They depend on the mass and the 
momentum of the outgoing lepton. Similarly, the func- 
tions Ro l( y 2 (ej q) are the nuclear response functions. The 
transition operators Cj(q), Lj(q), Ej(q), Mj(q) are the 
reduced Coulomb, longitudinal, transverse electric and 
transverse magnetic operators of angular momentum J. 
The response functions are calculated by inverting the 
Lorentz integral transforms 
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where a — <jr + iaj, and | ^i) ( 
the Schrodinger like equations 

(H-Eo-a) I %(<r,q)) = O i (q) | * ). 

The localized character of the ground state, and the 
imaginary part of a, give these equations an asymptotic 
boundary condition similar to a bound state. As a result, 
one can solve these equations with the hyperspherical 
harmonics (HH) expansion using the EIHH [ll|, EH G3 
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FIG. 1: (color online) LIT convergence for the leading multi- 
poles, as a function of the HH grand angular momenta K. 
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AV8' [3] 


AV18 


AV18+UIX 


AV18+UIX+MEC 


4 


2.09(-3) 


2.31(-3) 


1.63(-3) 


1.66(-3) 


6 


3.84(-2) 


4.30(-2) 


3.17(-2) 


3.20(-2) 


8 


2.25(-l) 


2.52(-l) 


1.91(-1) 


1.92(-1) 


10 


7.85(-l) 


8.81(-1) 


6.77(-l) 


6.82(-l) 


12 


2.05 


2.29 


1.79 


1.80 


14 


4.45 


4.53 


3.91 


3.93 



TABLE I: Temperature averaged neutral current inclusive 
inelastic cross-section per nucleon (in 10 -42 cm 2 ) as a function 
of neutrino temperature (in MeV). 



method. The matrix elements (^1^*2) are calculated 
with the Lanczos algorithm [2?3 |. 

In the supernova scenario one has to consider neutrinos 
with up to about 60 MeV. Usually, the leading contribu- 
tions in weak nuclear processes are the Gamow- Teller and 
the Fermi operators. Due to the total angular momen- 
tum and spin structure of the He nucleus, they are both 
strongly suppressed. In fact, the Gamow- Teller operator 
contributes only due to the small P- and -D-wave com- 
ponents of the ground state wave function. The same 
argument follows for the M\ operator. In addition, 4 He 
is an almost pure zero-isospin state [U HI] , hence the 
Fermi operator vanishes. Therefore, the leading contri- 
butions to the inelastic cross-section are due to the axial 
vector operators , , , L$ and the vector opera- 
tors CY,EY,LY (the latter are all proportional to each 
other due to the Siegert theorem). For the neutrino en- 
ergies considered here it is sufficient to retain contribu- 
tions up to 0(q 2 ) in the multipole expansion In Fig. 
[T] we present for these multipoles the convergence of the 
LIT as a function of the principle HH grand angular- 
momentum quantum number K . It can be seen that the 
EIHH method results in a rapid convergence of the LIT 
calculation to a sub-percentage accuracy level. Compar- 
ing with a previous work, [!|, we conclude that the 3NF 
does not affect much the convergence rate of these oper- 
ators. 

It is customary to assume that supernova neutrinos are 
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[MeV] 






(fe,e ) 


(^ e ,e + ) 


2 
4 
6 
8 
10 


1.47(-6) 
1.73(-3) 
3.34(-2) 
2.00(-l) 
7.09(-l) 


1.36(-6) 
1.59(-3) 
3.07(-2) 
1.84(-1) 
6.54(-l) 


7.40(-6) 
8.60(-3) 
1.63(-1) 
9.61(-1) 
3.36 


5.98(-6) 
6.84(-3) 
1.30(-1) 
7.68(-l) 
2.71 



TABLE II: Temperature averaged inclusive inelastic cross- 
section per nucleon as a function of temperature. 
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4.50(-2) 
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2.27(-l) 


1.80(-1) 
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9.26(-l) 
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4.56 


3.70 


8 


5.85 


5.43 


28.4 


22.9 


10 


21.7 


20.2 


103.8 


84.4 



TABLE III: Temperature averaged inclusive inelastic energy 
transfer cross-section per nucleon as a function of tempera- 
ture. 

in thermal equilibrium, so their spectra can be approxi- 
mated by the Fermi-Dirac distribution with a character- 
istic temperature T. In Table U we present the tempera- 
ture averaged total neutral current inelastic cross-section 
as a function of the neutrino temperature for the AV8', 
AV18, and the AV18+UIX nuclear Hamiltonians and for 
the AV18+UIX Hamiltonian adding the MEC. From the 
table it can be seen that the low-energy cross-section is 
rather sensitive to details of the nuclear force model (the 
effect of 3NF is about 30%). Sensitivity that gradually 
decreases with growing energy. In contrast the effect of 
MEC is rather small in our case, being on the percent- 
age level. Similar tendencies were also observed for the 
hep process • The small contribution of the MEC can 
be understood in the following way. The symmetry of 
the vertices in this low-energy approximation, dictate a 



symmetry between the two nucleons interacting via the 
meson exchange. The leading one-body multipoles have 
negative parity, as a result the MEC contributions to 
them is small. In comparison the MEC correction to 
the Gamow- Teller is of the same magnitude as the one- 
body current, however both terms become marginal with 
increasing momentum transfer. Although presented for 
the neutral current, these arguments hold true also for 
the charged currents since the response functions are re- 
lated by isospin rotation. 

In Table M and Table |TTT] we present (for 
AV18+UIX+MEC) the temperature averaged cross- 
section and energy transfer as a function of the neutrino 
temperature for the various processes. In both tables it 
can be seen that the charged current process is roughly 
a factor five more efficient than the neutral current 
process. Our results are of the same order of magnitude 
of previous estimates by Woosley et. al. [§[, though the 
differences can reach 25%. The current work predicts 
a stronger temperature dependence, with substantial 
increment at high temperatures. This indicates a 
different structure of the predicted resonances. 

Summarizing, we have given the first full microscopic 
study of v — a reactions, using a state of the art nu- 
clear Hamiltonian including MEC. The overall accuracy 
of our calculation is of the order of 5%. This error is 
mainly due to the strong sensitivity of the cross-section 
to the nuclear model, which could enhance the effect of 
uncertainties in the Hamiltonian. The numerical accu- 
racy of our calculations is of the order of 1%, and the 
cutoff dependence of the MEC is of the same order. With 
the present calculation, we make an important step in 
the path towards a more robust and reliable description 
of the neutrino heating of the pre-shock region in core- 
collapse supernovae, in which 4 He plays a decisive role. 
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